function [mu_I, L_I] = rnd_item_hyperparams(V)
[D,M] = size(V);

beta0 = 1;
nu0 = D;
mu0 = zeros(D,1);
W0 = eye(D);

V_ = mean(V,2);
S_ = fcov(V,2);

beta0_ = beta0 + M;
nu0_ = nu0 + M;
mu0_ = (beta0*mu0 + M*V_)/(beta0+M);
W0_ = inv(inv(W0) + M*S_ + (beta0*M/(beta0+M))*(mu0-V_)*(mu0-V_)');

W0_ = (W0_+W0_')/2;
L_I = wishrnd(W0_, nu0_);
Li_I = inv(beta0_*L_I);
mu_I = rnd_mvn(mu0_, Li_I);
